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Abstract. We consider the Derivative NLS equation with general 
quadratic nonlinearities. In [2] the first author has proved a sharp 
small data local well-posedness result in Sobolev spaces with a 
decay structure at infinity in dimension n = 2. Here we prove a 
similar result for large initial data in all dimensions n > 2. 



1. Introduction 

The general Cauchy problem for the semilinear Schrodinger equation 
has the form 

iu t - Au = P(u, u, Vm, Viz), t e R, x e R n 
u(x, 0) = u (x) 

where u : R n x R -> C m and P : C 2n+2 -> C m . Assuming that P is a 
polynomial containing terms of order at least k > 2 and higher, one is 
interested in studying the local well-posedness for this evolution in a 
suitable Sobolev space. 

In the simpler case when P does not depend on Vu, Vw this problem 
is rather well understood, and the Strichartz estimates for the linear 
Schrodinger equation play a proeminent role. 

Here we are interested in nonlinearities which contain derivatives. 
This problem was considered in full generality in the work of Kenig- 
Ponce-Vega [8]. Due to the need to regain one derivative in multilinear 
estimates, they use in an essential fashion the local smoothing estimates 
for the linear Schrodinger equation. Their results make it clear that 
one needs to differentiate two cases. 

If k > 3 then they prove local well-posedness for initial data in a 
Sobolev space H N with iV sufficiently large. However, if quadratic 
nonlinearities are present, i.e. k — 2, then the local well-posedness 
space also incorporates decay at infinity, namely 

H N ' N = {x N u e L 2 - D N u e L 2 } 
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with N large enough. The need for decay is motivated by work of 
Mizohata [10]. He proves that a necessary condition for the L 2 well- 
posedness of the problem: 



(2) 



iu t — Au = bi(x)Vu, t G 
u(x, 0) = u (x) 

is the uniform bound 



(3) sup 

xeR n ^es™- 1 ,R>o 



R 

Re I bi(x + roo) ■ uodr 
o 



< oo 



(4) 



The idea behind this condition is that Re b\ contributes to exponential 
growth of the solution along the Hamilton flow of the linear Schrodinger 
operator. 

If ac > 3 then naively one needs such bounds for expresions which 
are quadratic in u. If u G H N then u 2 G W N,X and the above integra- 
bility can be gained. However, if k = 2 then one would want similar 
bounds for linear expressions in u; but this cannot follow from square 
integrability, therefore one compensates by adding decay at infinity. 

In this work we consider the case of quadratic nonlinearities, precisely 
the problem 

iu t — Au = B((u, u), (Vu, Vm)), 
u(x, 0) = uq(x) 

Here B is a generic bilinear form which contains one differentiated and 
one undifferentiated factor, and also may have complex conjugates. 
Our results easily transfer via differentiation to the similar problem 
with two derivatives in the nonlinearity, 

iut -Au = B((Vu, Vu), (Vu, Vu)), 
u(x, 0) = u (x) 

The initial data u is assumed to be locally in Sobolev spaces H s but 
with some additional decay at infinity. In what follows A is a dyadic 
index. For a function u we consider a Littlewood-Paley decomposition 
in frequency 

u = J^u A , u x = S x u 

A>1 

where all the frequencies smaller than 1 are included in u\. Then 
following [2] we define the spaces VH S by 



(5) 



£ A2s IKI 
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where the dyadic norms T>L\ are defined in a manner somewhat similar 
to (3), namely 



This definition is consistent with the speed of propagation properties 
for the linear Schrodinger equation. Waves with frequency A have speed 
A therefore move about A within a unit time interval. Hence the linear 
Schrodinger equation is well-posed 1 in VH S . 

If one considers the low regularity well-posedness for (4), respectively 
(5) in VH S then a natural threshold is given by scaling, namely s c = 
| — 1 for (4) respectively s c — | for (5). However, it turns out that the 
obstruction identified in Mizohata's work is much stronger and leads to 
additional restrictions. In [3] Chihara obtains some better results on 
this problem, lowering the threshold for (5) to s = | + 4. However, this 
is still far from optimal. Indeed, a sharp small data result is obtained 
by the first author in a recent paper: 

Theorem 1 (Bejenaru [2]). If n = 2, the equations (4), respectively 
(5) are locally well-posed for initial data which is small in VH S for all 
s > s c + 1 . 

Although the above theorem was proved in dimension two, the same 
result can be derived in all dimensions n > 2. In addition, if some 
limited spherical symmetry is imposed on the data then the above 
exponents can be relaxed up to scaling, see [1]. 

The goal of this paper is to obtain a local well-posedness result for 
the same initial data space as in [2], but for large initial data. Our 
main result is 

Theorem 2. Let n > 2 and s > s c + 1. Then the equations (4), 
respectively (5) are locally well-posed for initial data in VH S . More 
precisely, there is C > so that for each initial data u £ VH S there 
is an unique solution 



In addition, the solution has a Lipschitz dependence on the initial data. 

We believe that this result is sharp for generic nonlinearities. How- 
ever, there are special cases when one is able to obtain stronger results, 

1 The operator norm of the evolution in T>H S will grow polynomially in time 
though. 




u E C(0,T;VH s ), 
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see for instance [4], [5], [6] ,[7], [12]. These results consider nonlineari- 
ties of the type B((u, u), Vtt) and B(Vu, Vw) in various dimensions and 
show that the regularity threshold for the initial data can be lowered 
all the way to the scaling. 

The step from small to large data is entirely nontrivial. The dif- 
ficulty is related to the infinite speed of propagation for the linear 
Schrodinger equation. Precisely, a large low frequency component of 
the solution produces an exponentially large perturbation in the high 
frequency flow even for an arbitrarily short time. Thus one needs to 
add the low frequency part of the data to the linear equation and only 
then do a perturbative analysis for the high frequencies. A somewhat 
similar analysis has been carried out before for related problems, see 
[8] , [3] ; however, in both works the full initial data becomes part of the 
nonlinearity, leading to considerable technical difficulties. What also 
differentiates the present work is that the perturbed linear equation is 
very close to Mizohata's necessary condition. 

Another interesting feature of this work is the choice of the function 
spaces for the perturbative analysis. It has been known for some time 
that the X s ' b spaces are not good enough in order to study the local 
theory for derivative NLS equation. This is due to the need to regain 
one derivative in the bilinear estimates, which leads to logarithmic 
losses. In [2] the first author introduced a refinement of the X s,b spaces 
which removes this difficulty In this article we provide an alternate 
modification of the X s ' b spaces which seems better suited for the study 
of variable coefficient equations. This is based on a wave packet type 
decomposition of solutions on the uncertainty principle scale. 

2. Scaling and the perturbative argument 

Differentiating once the equation (5) we obtain a equation of type 
(4). Thus in what follows we restrict our analysis to (4). We denote 

M = \\u \\ VH s 

Our equation is invariant with respect to the scaling 

(6) u £ (x,t) = eu(ex,e 2 t), u £ (x) = eu (ex) 

It is natural to seek to decrease the size of the initial data is by rescaling 
with a sufficiently small parameter e. However, there is a difficulty 
arising from the fact that we are using inhomogeneous Sobolev spaces. 
This is why the result in [2] does not yield a result for the large initial 
data. 
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We split the rescaled initial data u £ into low and high frequencies, 

U = U 0,<1 + U 0,>1 

It is not too difficult to estimate their size: 

Proposition 1. Assume that s > | . Then the components u £ Q<1 an d 
Uq yl of Uq satisfy the pointwise bounds 

(7) K,<ilU<» + IK>ilU~ ^ eM 
and the L 2 estimates 2 

(8) K,<ill2>tf<M, IIV^H^ <max{£,5 s -t}M, 
respectively 

(9) K>ill^ < eT^M. 

The low frequency component is large but has the redeeming feature 
that it does not change much on the unit time scale; hence we freeze 
it in time modulo small errors. On the other hand the high frequency 
component is small, therefore we can treat it perturbatively on a unit 
time interval. To write the equation for the function 

we decompose the quadratic nonlinearity B depending on whether the 
gradient factor is complex conjugate or not, 

B((u, u), (Vu, Vu)) = B ((u, u), Vu) + Si((u, u), Vw) 

Then v £ satisfies 

J iv t -Av-Av = B((v, v), (Vu, Vu)) + A t v + A^«,<i) 

(10) U(0) = < >1 

where 

Au = J B (K i < 1 ,ug,<i)>Vu) 
is a linear term which is included in the principal part, 

A lV = B((v, v), (Vm^j, Vu^)) + ^(K^x, m £ ,<i)> Vu) 
is a linear term which can be treated perturbatively and 

^K<l) = 5(K,<l,«0,<l)» (V<<1, V^<!)) + A^<! 

represents the time independent contribution of the low frequency part 
of the data. 

2 The constant in the second bound needs to be adjusted to e| lne| if s = ^ + 1. 
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To solve this we need some Banach spaces X s , T>X S for the solution 
v respectively Y s , VY S for the inhomogeneous term in the equation. 
These are defined in Section 3. 

Instead of working with the linear Schrodinger equation we need to 
consider lower order perturbations of it of the form 

(iv t -Av-a(t,x,D)v = f(x,t) 

I '•!.<•.()) !/(■>■) 

Given any M > 1 we introduce a larger class Cm of pseudodifferential 
operators so that for a G Cm we can solve (11). To motivate the 
following definition, we note that the imaginary part of the symbol a 
can produce exponential growth in (11). We want to be able to control 
the size of the growth in the phase space along the Hamilton flow of the 
linear Schrodinger equation, which leads to a (possibly large) bound on 
the integral of a along the flow. We also want the integral of a along 
the flow to give an accurate picture of the evolution. To insure this we 
impose a smallness condition on the integral of derivatives of a along 
the flow. This motivates the following 

Definition 1. Let M > 1. The symbol a(x,£) belongs to the class Cm 
if it satisfies the following conditions: 

(12) M = sup[ \a(t,x + 2t£,£)\dt<oo 

x,£ Jo 

(13) sup/" \d^d^a{t,x + 2t^,^)\dt<C a ^5 \a\ + \/3\>l 

x,£ Jo 

with 5e M < 1. 

Then our main linear result has the form: 

Theorem 3. Let s G R and a G Cm- Then the solution v of the 
equation (11) satisfies the estimate: 

(14) \\X{o,i]v\\vxs < e M (\\f\\-DY° + \\g\\vH>) 

To apply this theorem in our context we need to show that 
Proposition 2. For u G VH S as above set 

a(x,Q = £o(«,<i,«o,<i),0 

If e < e~ CM with C sufficiently large then a G C c „m with c n depending 
only on the dimension n. 
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In order to iteratively solve the equation (10) we want bounds for 
the right hand side terms. The main one is a bilinear estimate, 

Theorem 4. The following bilinear estimate holds: 

Tl 

(15) \\B((u,u), (Vv,Vv))\\t>ys < \\u\\ V x' \\v\\vx' s> 2 
For the linear term A 1 we use the estimates 

Proposition 3. Let A 1 be defined as above. Then 

(16) ||Aiw||y s < IKo^llU^lklU^ II^HdY* < II W 0,<1 ||^ ||x>JY s • 

and one for the time independent term: 
Proposition 4. If s > 1 then 

(17) WBiul^Vul^) + AuI^Wvhs <mm(e,e s ~%)M 2 

Using the results above we are able to conclude the proof of Theo- 
rem 2. We choose e = e~ CM with C large enough (depending on s). 
Then we solve the rescaled problem (10) on a unit time interval using 
the contraction principle in the space X s . 

We define the operator Ti by w = %_f to be the solution of the 
inhomogeneous Schrodinger equation with zero initial data: 

iw t — Aw — A(x, D)w = f 
w(x, 0) = 

We also denote by T 2 g the solution to the homogeneous equation: 

iw t — Aw — A(x, D)w = 
w(x,0) = g 

With these notations the equation (10) can be rewritten in the form 
v = Tv, Tv = T 2 u £ Q>1 + %.(B(v, Vu) + A x v + B(w £ , Vw £ ) + Aw £ ) 
We define the set 

K = {w E VX S : \\w\\vx> < < a < s - - 

and prove that T : K — > K and that T is a contraction on K. This 
give us the existence of a fixed point for T which is the solution of our 
problem in the interval [0, 1]. 

To prove the invariance of K under the action of T we use the results 
in Theorem 3: 

\\Tv\\ VX s < e^fllu^HDff. + \\B(v,Vv)\\ VY s + 
+ ll s «<iV«o,<i) + Auo,<i\\vh') 



(18) 



(19) 
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Using the estimates (9), (15), (16) and (17) this yields 

\\Tv\\vx> <M 2 e c - M (e s -i + \\v\\ 2 VXa ) 

If v G K then we use the smallness of e to obtain 

\\Tv\\ VX s < e a M 2 e CnM (€ s -^- a + e a ) < e a 

which shows that Tv G K. 

To prove that T is a contraction we write 

Tui - Tv 2 = Ti(S(ui, V(ui - ua))) + Tx{B{v x - v 2 , Vu 2 )) 

and estimate in a similar manner 

||Tui - Tw 2 ||x>x« <(||ui||x>a-« + IM|x>xOIK - V 2 \\vX* 

<2e CT ||x[o,i]( w i -U2)||da-« < ^||X[o,i]( u i ~ V2)\\vx° 

3. The function spaces 

Let (xq,£o) G M. 2n . To describe functions which are localized in the 
phase space on the unit scale near (xo,£o) we use ^ ne norm: 

H»% := {/:<£>- 0*7 G L 2 , <* - G L 2 } 

We work with the lattice Z ra both in the physical and Fourier space. 
We consider a partition of unity in the physical space, 

0*o = 1 0xo = <K X ~ x o) 

where is a smooth bump function with compact support. We use a 
similar partition of unity on the Fourier side: 

Let if be a Hilbert space. Let V 2 -ff be the space of continuous H 
valued functions on R with bounded 2-variation: 



|«" 2 



2 V 2 H = sup J2\\u(t i+1 ) -u(U)\\ 2 H 

where T is the set of finite increasing sequences in R. We abuse nota- 
tions and consider V 2 to be the closure of the space of smooth functions 
in this norm, namely we use V 2 instead of V 2 fl C. The V 2 spaces are 
close to the homogeneous Sobolev space H 2 in the sense that 

(20) C V 2 C fl* 
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Let U 2 H be the atomic space defined by the atoms: 

i i 

for some (U) G T. We have the inclusion U 2 H C V 2 H but in effect 
these spaces are very close. There is also a duality relation between 
V 2 H and U 2 H, namely 

(21) (DU 2 H)* = V 2 H 

We can associate similar spaces associated to the Schrodinger flow 
by pulling back functions to time along the flow, e.g. 

II II II it A ii 

INIv^L 2 — \\ e U\\v^L 2 

This turns out to be a good replacement for the X '^ space associated 
to the Schrodinger equations. Such spaces originate in unpublished 
work of the second author on the wave-map equation, and have been 
succesfully used in various contexts so far, see [9] and [11]. 

In the present paper we consider a wave packet type refinement of 
this structure. We begin with spatial localization, and introduce the 
space X of functions in [0, 1] x W 1 with norm 

\Hx= ll^o^Nl^i 

x ez 2 

It is easy to see that this is a stronger norm than the V^L 2 norm, 

(22) \W\\vli? ^ \\u\\ 2 x 

The X norm is usually applied to functions which are also frequency 
localized on the unit scale. This is consistent with the spatial localiza- 
tion. Precisely, we have the straightforward bound 

(23) \\M D Mx < \\u\\x 

For a function u : [0, 1] x R n — > C we decompose 

u = Yl u & = ( Pto( D ) u 

and set 

II*. = E<£o> 2 1K,llx 



An immediate consequence of (22) is that 
(24) hWvgL* ^ h\\x° 
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We also denote by X\ the subspace of functions in X° which are 
localized at frequency A. It is easy to see that X s has an I 2 dyadic 
structure, 

\\uf xs « ^A 2s |K||^ A 

A>1 

This is the most compact definition of X s . Some equivalent formula- 
tions of the norm on X turn out to be more helpful in some estimates. 
We dedicate the next few paragraphs to such formulations. The first 
such formulation simply adds regularity and decay: 

Proposition 5. Let N e N. Then 

(25) IKollx- H^oe itA u&llv 2H "." 

X06Z* X °' (0 

Proof. Denote v^ = e ltA u^ . Then it suffices to show that 

z2 HxoV^WL h n,n < V ||0 

xn, En 

xoez* xoez* 
This in turn follows by summation from 

110*0^0 II van"-? ~ y2i x o-yo)~ N \\<i>voVf \\v*L* 

y ei n 

We can translate both in space and in frequency and reduce the problem 
to the case when x = and £ — 0. Then we need to show that 

11^00^0 1| y2L2 + \\D N {4>ov )\\ V 2 L 2 < ^2 (yo)' N \\ ( Pyo v o\\v 2 L2 



The derivatives which fall on vq can be truncated at frequencies larger 
than 1. Then a more general formulation of the above bound is 

(26) ||^oXo(-D)uo||v^ < (y°)~ N II ^yo v o II v 2 L2 

where both ip and Xo are bump functions concentrated at 0. We write 
*PoXo(D)v = Y ^oXo(D)(f)y v 



It remains to show that 

||Voxo(£>)<M^->l2 <{yo)~ N 

which is straightforward since Xo(D) has a bounded and rapidly de- 
creasing kernel. □ 
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The next equivalent definition of our function spaces relates the 
Schrodinder evolution to the associated Hamilton flow, 

(zo,fo) -> (x - 2if , fo) 
We begin by the linearizing the symbol of —A near £o, 

e = l 5o (o + o((e - Co) 2 ), ^ (o = e 2 + 2^0 

The evolution generated by L is simply the transport along the Hamil- 
ton flow, 

e itL ^ u (x) = e it $u(x-2t£ ) 
Our next characterization of the frequency localized X norm asserts 
that we can replace —A by L^ . 

Proposition 6. We have 

(27) ll^o llx~ E W^ e ~ itLsou ^W 2 v^ 

x a ez 2 

and 

(28) IKIIx~ E ll^oe- itL£o «&II^H^ 

Proof. Denoting as before = e ltA u^ we can write 

e-^o^ o = e ^ D2 - L to\ = X^ (t,D)v (0 

where the symbol is a un ff bump function in £ around £o and 

smooth in t. There is also a similar formula with L^ and D 2 inter- 
changed. Hence for (27) it suffices to show that 

Without any restriction in generality we can take £o = 0. Using a 
Fourier series in t we can also replace Xft(A0 by an expression of the 
form a(t)x£ (Q with a smooth. It remains to show that 

E WKX^( D Xo\\v 2 L 2 ^ ^2 H xo v £,o \ \v 2 L 2 
x a ez 2 x a ez 2 

But this follows as in Proposition 5 from (26). 

Finally, the proof of (28) uses the same argument as in Proposi- 
tion 25. □ 

For each £ we define the set of tubes Tg of the form 

Q = {(t,x) : \x - (x Q - 2t£ )\ < 1} Xq E Z n 
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They are the image of unit cubes centered at xq at time along the 
Z/£ flow. These tubes generate the decomposition: 

[o, i] x r = |J q 

Now we can state our last equivalent formulation of the X norm for 
functions localized at frequency £o i n terms of a wave packet decompo- 
sition associated to tubes Q G Tg : 

Proposition 7. Let u^ G X . Then it can be represented as the sum a 
rapidly convergent series 

(29) u^t, x) = J2 e^e*® E <j(t)^(x - x Q - 2^o) 
with x 3 uniformly bounded in H N,N , supported in 5(0,2) and 

E W a Q\\v 2 %i~ N \\ u to\\x 

Proof. The pull back to time of the above representation using the 
L^ flow is 

e-^o M& (t,x) = ^e^' ^(t) X j (x-x ) 

We denote t> = e~ lx ^°e~ ltL ^u^ and use (28) for the X% norm. Then 
it suffices to show that for fixed Xq we can represent 

<t>x VQ = E a3 ( x ~ x °"> 

3 

where x 3 are uniformly bounded in H N,N and 

\\ a3 \\v* ^ II 0Z(,^o|| y2 H CN,CN 

Xq 

Without any restriction in generality we take xq = 0. We take x 3 to 
be the Hermite functions with the H N ' N normalization (see e.g. [13]). 
Then the V 2 functions a 3 are the Fourier coefficients of 4> Xo vq. They 
decay rapidly due to the additional regularity of (p xo v . 

Finally, to insure that the x/ s have compact support we can truncate 
the Fourier series outside the support of <p xo . 

□ 

We follow a similar path to define the Y s structure: 
II*. = E^) 2 1Klly 
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where Y is defined by: 

11*= E WKe- U ° 2 f\\l u!Ll . 

All the equivalent definitions for the X have a counterpart for Y by 
simply replacing the V t 2 structure by DUj 1 one. 

There are two key relations between the X s and the Y s spaces. The 
first one is concerned with solvability for the linear Schrodinger equa- 
tion: 

Proposition 8. The solution u to the linear Schrodinger equation 
iu t - Au = /, u(0) = u 

satisfies 

(30) \\u\\x> ^ ||«o||h« + WflW* 

This is stated here only for the sake of completeness, as in the next 
section we prove the stronger estimate in Theorem 3. 
The second relation is a duality relation: 

Proposition 9. We have the duality relation 

(Y s y = x- s . 

Proof, a) We first verify that X~ s C (Y s )*. For this we need the bound 

.i f 

ufdxdt 



We decompose 



/ / ufdxdt = / / u^f^dxdt 
Jo J £ gZ „ Jo J 



Due to the definition of the X s and Y s norms, it suffices to show that 



o 



u^f^dxdt 



< H\\x\\h\\y- 



For this we write 

-1 r rl 



ufdxdt = J eit ^H e itA f^dxdt 

= E / / K^ tA H ^T^dxdt 

xo ez* Jo J 
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and use the duality relation (21) together with the definition of the X 
and Y norms. 

b) We now show that (Y s )* C X~ s . Let T be a bounded linear 
functional. Then we have 



\Tf\ < ||/||y. < (E^IIAII 



< 



E J2^ 2s \\Ke UA h\\ 2 DU^ 



By the Hahn-Banach theorem we can extend T to a bounded linear 
functional on the space l 2 ^ s DU 2 L 2 . By (21) this implies that we can 
represent T in the form 

T f = E E f I v X0 ,zKe UA Udxdt 

where 

V s /t\-2s||„. 1 1 2 <- n/71112 



E E^iKdi^im 12 



Due to the above representation we can identify T with the function 

Ut= E E ^( £) ) e_<tA ^o^o,C 

It remains to show that 

Kllx- < E E^^ll^o.dl 2 /^ 2 

This reduces to the fixed £ bound 

II E M D ) e ~ ttA( />xo v xodx £ E IK>,dlv 2 L 2 

with a modified </?£. Using the definition of the X norm we rewrite this 

as 

E ii E ^w^cWso^wiiv^ ^ E n^o.di 2 /^ 2 - 

But this follows by Cauchy-Schwartz from the rapid decay 
\\<f)y ^{D)(f) X0 \\ L 2^ L 2 < (x - y )~ N - 

□ 
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We still need to add the decay structure to the X s , respectively the 
Y s spaces. Given A > 1, we roughly want to ask for I 1 summability of 
frequency A norms along collinear cubes of size A. We define the VX X 
norm by 



oo 



(31) 11^11^= sup sup S~] ||x(A 1 (x - x ) - ku)u\\ Xx 

where x is a compactly supported bump function. We note that this 
norm can only be meaningfully used for functions at frequency A. Sum- 
ming up with respect to A we also set 

ii ii 2 — V s ii ii 2 

A>1 

In a completely similar way we can define T>Y\ and T>Y S . 

A useful tool in our analysis is a family of embeddings which corre- 
spond to the Strichartz estimates for the Schrodinger equation. 

Proposition 10. Let p and q be indices which satisfy 
2 n n 

- + - = -, 2 < p < oo, 2 < g < oo 
p q 2 

Then we have the embedding 

X° C LfLl 



By (24) these embeddings are a direct consequence of the Strichartz 
estimates for the Schrodinger equation, see [9] , the proof of Proposition 
6.2. By Sobolev embeddings we also obtain bounds with larger p, q for 
frequency localized solutions. 

Corollary 1. Let p and q be indices which satisfy 

2 Tl Tl 

- + -<-, 2<p<oo, 2<g<oo 
p q 2 

Then we have the embedding 

n Tl 2 n „ 

S <fM X° C p--I%Ll 

where S <lx can be replaced with a multiplier localizing in frequency to 
an arbitrary cube of size [i. 

In this article we use only the case p = q, more precisely 

(32) I°CL^, S^X C ^L n+2 
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Finally, we introduce modulation localization operators, which we 
can define in two equivalent ways. First is as multipliers, 

M <a u = s <a (r-£ 2 )u. 

The second is obtained by conjugation with respect to the Schrodinger 
flow, 

e aD \M <a u)(t) = s <a (D t )(e itD2 u(t)) 

The operators s <cr (D t ) are bounded on V 2 . On the other hand for 
the remainder, by (20), we have a good L 2 bound: 

\\s >a (D t )a\\ L 2 < (7~5|| a ||y2. 

By the second definition of the modulation localization operators above 
we obtain 

Proposition 11. a) The operators M <a are bounded on X° . 
b) The following estimate holds: 

\\M >a u\\ L 2 < a~^\\u\\ X o 



We note that in effect, by using both inclusions in (20), one can 
relate our X s spaces to the traditional X s ' b spaces, namely 

(33) X 3 '^ 1 C X s C X s '5'°° 

4. Linear estimates 

This section is devoted to the study of the linear equation (11) with 
a G Cm- Precisely, we aim to prove Theorem 3. We denote 

L a = idt — A — a(t, x, D) 

We begin our analysis with a heuristic computation. Suppose we 
have a solution u x °^° to 

(id t - A-a(t,x,D))u = 

which is localized on the unit scale near the bicharacteristic t — > (xq + 
2t£ , Co)- Then we can freeze the symbol of a(t, x, D) along the ray and 
write 

a(t, x, D)u = a(t, xo + 2£ , £,o)u + error 
Thus it 10 ''' approximately solves 

(id t - A - a(t, x + 2£ , &))«*,,,& ~ 
This implies that we can represent u x °^° in the form 

u Xo '^°(t) ~ a( - s ' Xo+2s ^ Q '^ ds v x "^° 
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where v x °^° solves the equation 

(id t - A)v xo4 ° = 

Hence along each wave packet we can use the above exponential to 
approximately conjugate the variable coefficient equation to the flat 
flow. 

Using this idea we produce wave packet approximate solutions for 
the equation (11). By orthogonality these combine into general approx- 
imate solutions to (11). The exact solutions are obtained via a Picard 
iteration. 

We first consider the regularity of the exponential weight. By (12) 
we have 



Jo 



\a(t,x + 2t^o,^ )\dt < M 



which implies a W ' bound for the exponential, 



(34) 



^_ e J* a(t,x +2t^ 4o)dt 

dt 



< e 



M 



Since W 1 ' 1 C V 2 this yields a similar V 2 bound, 



(35) 



J*a(t,x +2t€o,£o)dt 



V 2 



< e 



M 



We continue our analysis with the localized equation 

(36) L a u = f X0 £ , u(0)=g XOtCo 

According to the above heuristics, an approximate solution for this 
should be given by Duhamel's formula, 

M ^o,5o _ e f£a(s,x +2s£o,€ods) e -itAg^ 
ft 

+ / e ^ a(T ' Xo+2r?0 ' 5o)dr e^^ s ) A / ;c0i5o (s)rfs 
Jo 

We prove that this is indeed the case: 

Proposition 12. The function u x °'^° is an approximate solution for 

(36) in the sense that 

(37) ||e*< A ^|| y2 ^ < e M (\\g Xo ^\\ H , N + W^U^du^") 
and 

|| e «A ( LaU *o*> _ fxQ ^ || < 

(38) x ° 

3e M ( \\g XQ £ || H N+2n+l,N+2n+l + 1 1 e*' A f XQ £ || D y 2 ff JV+2n+l,JV+2n+l ) 
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Proof. For the first bound we shorten the notation 

a (s) := a(s, x + 2s£ , £o) 

and compute 

Jo 

The first factor is estimated directly by (35). Setting 

F(t)= fe^&Mds 
Jo 

and integrating by parts we write the second term in the form 

F (t) ~ f ^-ef> it>T)dT F(s)ds 
Jo ds 



We have 



\ F \\v 2 H N '" ~ \\ eiSA fxo&X S ) 



■ ■ . i>r'n s ) 



while the V norm of the second part is controlled by its W ' norm, 



namely 

" d 



f A e />(or)dr F(s)ds 

Jo ds 



y2 H N,N 



< 



4 f ^ e l>0(r) d r F{s)ds 

dt J ds 



L1H N,N 



< 



M*)l 1 



Jo 



A e /.'«o(T)dr 

ds 



ds dt\\F\\ LCOH N, 



< 



KM 



(1+ / \ao(s)\e^^^ds)dt\\F\\ LooH 



< f 1 \a (t)\efoM^dt\\F\\ LooHN , 
Jo x ° 



N 

■Co 



< v M II z?ll 



< v M \\ Z?ll 
- e IF II V 2 H N ' N 



This concludes the proof of (37). 

It remains to prove (38). A direct computation yields 



(i dt _ A - a(t, x, D))u x ^ - / X0)& = &(f , D)u x °^° 



where 



&(*, ^> = + 2^o, Co) - a(t, C) 

By (37) it suffices to show that 

\\e itA b(t,x,D)u X0 't°\\ LlH N,N < 5\\e itA u X0 't°\\ LooH N + 2 n+ i,N + 2n + i 
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Since the flat Schrodinger flow has the mapping property 
\\e itA f\\ H > 



jN.N ~ \\J\\ H N,N 

this is equivalent to 

(39) \\b(t,x,D)u x ^\\ LlHN ,„ <6\\u**\\ L „ H »w+* 
We begin with a straightforward consequence of the S 00 calculus: 

Lemma 1. Let k be a nonnegative integer and c be a symbol which 
satisfies 

\d^c(x, 1 < c a , p ((x - so, £ - &)>*, \a\ + \/3\>0 
Then for all N G N we have 

\\c(X } D)u\\ H N,N < \\u\\ H N+k,N+k 

In order to use this lemma for the operator b above we need the 
following Sobolev embedding: 

Lemma 2. Let R > 1 and c G W 2n ^{B R {xQ^ Q )). Then 
\c(x,0\< J2 H^IUw^o)) 

0<|a|<2n 

respectively 

\c(x,£)-c(x ,So)\< W^hHB^xo^o)) 

l<|a|<2n 

As a consequence of this we obtain 
Lemma 3. Let c G W^' 1 (R 2n ) . Then 

\c(x,Q\ f v l^cfofll 



respectively 



\c(x,g)-c(x ,Zo)\ ^ f ^ 1^(^01 

l<H<2ra 



< > ; 77- 1 T' ; l,m 



Applying these inequalities to the symbol b above we obtain point- 
wise bounds for b 

1^,01 t E5-,ia°°(^.0l 



((x - - 2t&, { - { ))2»+i (i - 1 - 2t&, 5 - {„)) 2 »+! 
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and also for its derivatives, 

Then by Lemma 1 it follows that 6 satisfies the bound 

with iVo sufficiently large. Integrating in t and changing coordinates in 
the integral this gives 

\\b(t,X,D)\\ LXH N+2n+l,N+2n+l_^ LlH N,N^ ^ < 

• 1 Wo 



U ^ ((i-io.{-«o)> 2 " +1 

By (13) we can bound each time integral by 5 and the remaining weight 
is integrable in x and £. Hence (39) follows. 

□ 

Next we produce approximate solutions for the frequency localized 
data, 

(40) L a u = fz , u(0)=g (o 

We denote the approximate solution by u^°; the notation u^ continues 
to be reserved for a frequency localized part of a function u. 

Proposition 13. There is an approximate solution to the equa- 
tion (40), localized at frequency £ ; with u(0) = g^ , which satisfies the 
bounds 

(41) ll^°IU<e M (|| fe |U 2 + ||/, || y ) 
respectively 

(42) ||S € (L a ti*» - f 6o )\W < S(d ~ Zo)- N e M (\\g 6o \\Li + ||/ & ||y) 
Proof. We decompose 

9(o = ^ ] 9xq,£o f%o = ^ ] fxp,£o 
xo xo 

and solve the problem (36) for which we have the estimates (37) and 
(38). We define our approximate solution to be the sum of the approx- 
imate solutions u x °'^°: 



Xo 
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Using (37) we obtain: 



< 



yo 

yo V x J 



\ XQ,Zo 

yo \ xo j 

< llp- itA 7l X 0'^H 2 
Xo 

<p M Vlln II 2 i L-flAf ||2 

~ C / ^ Ili/XO^O ||jjJV+2n+l,JV+2n+l "T ||C J X ,£ \\ D{ j2 H N+2n+l,N+2n+l 

<e M (lkoll! 2 + llAoll^ ) 

We continue now with the estimates for the error, using (38) instead 
of (37). We have 

\\Ss(L a ut° - fail = E Uye itA Ss(L a uZ° - h )\\lu^ 



v 



< 



E E UyS^ tA (L aU ^ - f xo ,a\w 



y \ xo 



2 



<(v-*o-2teo, 

<<£ - e ) n - 2JV E 11^^° - /xo^b)ll^ 

Then (42) follows from (38). □ 
The next stage is to consider data which is localized at frequency A, 
(43) L a u = f x , u(0) = g x 

Summing up the frequency localized solutions we obtain as above a 
dyadic approximate solution: 

Proposition 14. There is an approximate solution u x to the equation 

(43) , localized at frequency \, with u(0) = gx, which satisfies the bounds 

(44) \\u X \\x x <e M (\\gx\\ L2 + \\h\\yJ 
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respectively 

(45) \\S,(L a u x - f x )\\^ < (mini^^y 6e M (\\g x \\ L2 + \\f x \\ Yj ) 

The construction of the functions u x involves only the constant co- 
efficient Schrodinger flow at frequency A. This has spatial speed of 
propagation A, therefore it can spread by at most 0(A) in a unit time 
interval. Thus the above construction can be trivially localized on the 
A spatial scale, leading to the bounds 

(46) h%x x <e M (\\g x \\ VLl + \\f x \\vY x ) 
respectively 

/ \ \ N 

(47) \\S,(L a u'-f x )\\ VY ^<(mm{^-}) 5e M (\\g x \\ VLl + \\f x \\ VYx ) 

After an addition dyadic summation we obtain a global parametrix: 

Proposition 15. There is an approximate solution u to the equation 
(11), withu(0) = g, which satisfies the bounds 

(48) \\u\\vx°<e M (\\g\\ VHS + \\f\\ VY s) 
respectively 

(49) \\(id t - A - a(t,x,D))u - f)\\vy. < 5e M (\\g\\ VH s + \\f\\ VY s) 

Of course the similar result without decay is also valid. 

If S e _M then the constant in (49) is less than 1. Then one can 
iterate to obtain an exact solution u to (11) which still satisfies (46). 
Theorem 3 follows. 

5. LOW FREQUENCY BOUNDS 

Here we prove Propositions 1,2,3,4. 

Proof of Proposition 1. For simplicity we denote / = u . The point- 
wise bounds (7) follow from Sobolev type estimates and scaling, 

71 

ll/lL-=e||/|U-<e||/|| OT ., *> 2 

We now move to the L 2 bounds. The effect of scaling on the VL\ 
norms is easy to compute, 

(so) \\nvLi = ^-nf\\ VLlx 

Hence the first bound in (8) can be rewritten in the form 
(51) e^ll/^-iH^^H/llM. 
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We have 

/<£-! — E fx 

1<A<£-! 

For each such A we have three relevant spatial scales, 

e < A -1 < A 

The middle one arises due to the uncertainty principle; namely, this is 
the scale on which fx is smooth. Then we can write the sequence of 
inequalities 

\\h\\vi? e ^ {eXYh^WfxWvL^ 

= e^X^WfxWvLl 

All these bounds are obtained by comparing tubes with the same ori- 
entation uj. The first step uses Holder's inequality to switch from the 
e scale to the A" 1 scale; the VL°° norm is defined in the same way 
as the DL 2 norm. The second takes advantage of the localization at 
frequency A. Finally the third step uses Holder's inequality to switch 
from the A -1 scale to the A scale. Summing up with respect to A we 
obtain (51), 

ll/<e-i||MS< E Wh\\vLl<e^ £ \Hh\\vLl<S^\\f\\ VH s 

for s > n/2. 

Consider now the second part of (8). The analogue of (50) is 
\\^r\\vLl=e 2 -nVf\\ VLlx 
therefore we have to show that 

e 2_5 ||V/< e -i|| M a <max{e,£ s -t}||/||^ s 
We proceed as above, 

l|V/<e-i||^< E M\h\\vLl<^ 1 E Af+1 H/A||^ 

1<A<£-! l^A^E- 1 

The bound for the last sum is straightforward and depends on the 
relative positions of s and | + 1. 
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It remains to consider the high frequency estimate (9), for which we 
only need the scaling relation (50) and Holder's inequality: 



ll/>llll>if s = ^2 X2S \\(f £ )^\\T)Ll = £2 n ^2 X ' 2S \\fe- 1 x\\iLi 

\>1 \>1 



A>1 

□ 

Proof of Proposition 2. It suffices to consider the a(t,x,l;) = «o <i£ 
part of the symbol, since Vmq^i is obtained from a by differentiation. 
The main ingredient of the proof is 

Lemma 4. Let s > | . Then the following estimate holds: 

/oo 
\f(x + tuo)\dt < \\f\\vH° 
-oo 

Proof of Lemma 4- For each dyadic component fx of / we have: 

/oo _ /■A(fe+1) 

\f x (x + tu )\dt = J2 I/aK + tcu )\dt 

-oo k& Jxk 

( r x(k+i) 
/ \ft(x + tuj )\dt 



< At 



A2 ll^flA-H^-^-fcw^^/AlU 2 



£ A ? ||/a||©L? 



A 



At the second step we have used Holder's inequality and at the third 
we have used the frequency localization. The summation with respect 
to A gives the desired result in (52). 

□ 

We return to the proof of the proposition. A change of variables 
combined with Lemma 4 and the first part of (8) gives 

f 1 _ f m 
sup / \a(t,x+2t£,£)\dt = sup / {u^^x+tuo^dt <||mo,<iII^ s ~ M 

x,£ JO JO 

where oj — 
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For x derivatives of a we use the second part of (8) instead: 



25 



JO 



\d x a(x + 2t£,0\dt< \\Vu £ 0<1 \\ VH s <m&x{e,e s -%}M 



For £ derivatives we use the pointwise bound (7): 

\dta(x + 2t£,Z)\dt= [ \u e ^ 1 {x + 2tC)\dt<eM 
'o Jo 

Since u e Q K1 is supported at frequencies < 1, we obtain similar bounds 
for higher order derivatives. □ 



Jo 



Proof of Proposition 3. We neglect the gradients applied to u £ <v Also 
for the estimate involving v we retain only the stronger bound 3 involv- 
ing Vv. Then we need to show that 

( 53 ) \\ u o,<i v \\y° < ll«o,<ilU 00 ll u ll-x-« 
respectively 

(54) ll«o,<iVu||y« < K,<i||l~|MU 3 - 

In (53) we use orthogonality with respect to unit frequency cubes to 
reduce it to 

IK<i^olk ^ K,<ilU~IKIU 
By duality this becomes 



^ IK<ilU~lkolUIKJ x 



u £ Q<l v^ w^ dxdt 
Hence it suffices to show that 

\v^ \\w^\dxdt < IMUIKJx 



/ 



By orthogonality with respect to T^ tubes (see Proposition 7) it re- 
mains to verify that given a £ tube Q we can integrate a bump function 
on Q, 



I 



l Q dxdt < 1 



which is trivial. 

Similarly, (54) reduces to the frequency localized bound 



ul^v^w^dxdt 



^ 161 K,<ilk~IMIHKolk 



3 This is because we do not differentiate between frequencies 1 and less than 1 in 
what follows 
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We use a modulation decomposition of v^ and u>_g at modulation 
£ 2 /8. Due to the Fourier localization, the integral corresponding to the 
low modulation parts vanishes, 



u o,<i M <eg /4V '& M^/tW-frdxdt = 
On the other hand we use Proposition 11 to estimate 



/ 



^ 1 1 «0,< 1 1 U°° 1 1 M >$/4 V Zo I U 2 1 1 1 1 W ~to I U 2 
^ I Co | _1 1 K,< 1 1 1 L~ 1 1 V£ \\ X o\\W-z \\ x0 



□ 



Proof of Proposition 4- This is a direct consequence of the estimates in 
(8) and (7). □ 

6. The bilinear estimate 

Here we prove Theorem 4. After a Littlewood-Paley decomposition 
it suffices to prove a bound for the high-low frequency interactions 

(55) \\u^v x \\ VYx < A~V^( m ^IMIx>xJM|x>x A 1 < n < A 
respectively 4 for high-high frequency interactions 

(56) WS^UxVx)^ < A n ~V f IKI|wdMbx A 1</"<A 

and the similar bounds where one or both of the factors are replaced 
by their complex conjugates. 

Both bounds can be localized on the 2 J spatial scale. Thus the two 
D/s can be factored out in the first bound, and neglected in the second 
one. Then, by the duality result in Proposition 9, (55) can be rewritten 
as 



(57) 



J u^vxwxdxdt <A 1 /i2(ln/i)2 



2 "u^\\vxJ\vx\\x x \\wx\\x x 



On the other hand, in (56) we can replace the I 1 summation by an 
l 2 summation on the \x spatial scale by losing a \2ji~2 factor. Thus it 
remains to show that 

\\S^{uxVx)\\y^ < A n_ ^ _ ^|KIU A IMU A 

By the duality result in Proposition 9 this is equivalent to 

uxvxw^dxdt < \ n ~^' L ~ Tl \\ux\\x x \\vx\\x x \\w^\\x l , 



I 



Strictly speaking we should consider products of the form u Xl vx 2 with Ai 
but this makes no difference 
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which is easily seen to be weaker than (57). 

It remains to prove (57) where we allow 1 < fi < A in order to 
include both cases above, and where we allow any combination of com- 
plex conjugates. The seemingly large number of cases is reduced by 
observing that the bound rests unchanged if we conjugate the entire 
product. Hence we can assume without any restriction in generality 
that at most one factor is conjugated. Hence it suffices to consider the 
following three cases: 

(i) The product u^v \W\ with 1 < /j, < A. This is the main case, where 
all three factors can simultaneously concentrate in frequency near the 
parabola. 

(ii) The product u^v\W\ with 1 < /i C A. Because the high fre- 
quency factors cannot simultaneously concentrate on the parabola, the 
estimate turns essentially into a bilinear L 2 estimate. 

(iii) The product u^vxwx with 1 < /i < A. This is very similar to the 
second case. 

Case 1: Here we prove (57) exactly as stated, for 1 < fi < A. This 
follows by summation with respect to £, i] in the following result: 

Proposition 16. a) Let 1 < < A and £,r] G R™ with 
|£|~A, |£ + ??|~A. 

Then we have 



(58) 



S v uS^vS^ +ri wdxdt 



~ /^(A + min{|£ ■ rj\, |f Arj\})^ 



b) Assume in addition that /i < A 2 . Then we have 



(59) 



J S v uS^vS^ +v wdxdt 



< \\S v u\\v lX JS ( v\\x ( \\S i+v w\\x s+v 

A5/i"5(A+ ^-7/1)5 



We first show how to use the proposition to conclude the proof of 
(57). We decompose u^, v\ and ui\ in unit frequency cubes, 

M«*» l€l»A ICI-a 

and use the corresponding decomposition of the integral in (57). 
We consider two cases. If ji > A2 then we use (58) to estimate 



UfjV\W\dxdt 



|€|«£m«#. ^(A + min{|e^|,|eA77|})5 
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By the Cauchy-Schwartz inequality this is bounded by 



1] ll-VJU-W^II^+J l 

and further by 



\uJxJvx\\x x \\w x \\xA sup V* . r|< , ( , ... 

1 ici«a ,v M A + mm {l£-% K a 77I}) 



which gives (57) since 

1 u n m u 

' _ 

li(\ + mm{\Z-r)\,\£Ar)\}) " /i 2 A 



E 



and /i 2 > A. We note that the bound improves as /x increases. 

If /i < A2 then the argument is similar but using (59) instead of (58). 
This concludes the proof of (57). 

Proof of Proposition 16. We decompose each of the factors in wave 
packets, 

S v u = ^2 U P' S Z V = ^2 v ®i s z+v w = S Wr 

PeT v QeT 5 ReT 5+v 

We first prove a bound with Q, R fixed. 

Lemma 5. For £ and 77 as above the following estimate holds: 



(60) 



^ u P v Q w R dx 
Per v 



Proof. By Proposition 7 we can assume without any restriction in gen- 
erality that 

Mp = a P (t)e-^ 2 xp(a; - 2t77), t; Q = ag^e"^^* - 2i£), 

«J5 = a fl (*)e i *«- H ') 2 Xfl(a:- 2^ + 77)) 

with the x's being unit bumps and the a's in V" 2 . Then the integral 
has the form 

e 4it ^a Q {t)a R {t)J2ap(t) [ xp{x-2t V )x Q {x-2tOxR{x-2t^+ V ))dxdt 
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The tubes Q and R differ in speed by 77, therefore they intersect in 
a time interval / of lenght at most /j -1 . The tubes P and R differ in 
speed by £, therefore they intersect in a time interval of lenght at most 
A -1 . Thus there are about A/j^ 1 tubes P which intersect both Q and 
R. 

For fixed P the x integral above is a smooth bump function on a A -1 
interval. Then we can express the above integral in the form 



e 4U ^a Q (t)a R (t)Y,ap(t)b P (t)dt 



where bp are smooth bump functions on essentially disjoint A 1 inter- 
vals. For the sum with respect to P we can estimate 



^(OMOIIi^A-^llapl 



and 



$>p(i)MOIIy 2 <£ll 



We consider two possibilities. If |£ • 77 1 < A then we use Holder's in- 
equality to bound the integral by 



\dp 1 1 L°° 1 1 1 1 L°° 1 1 a R 1 1 L°° 



If |£ ' ^1 ~ ^ then we use the algebra property for V 2 . It remains to 
prove that 



a(t)e ita dt 
After rescaling this becomes 



< 



fi 2 \a\ 2 \\cl\\v 2 



4$V 



a(t)e lta dt 



^ W\ 2 ll a llv 2 



This follows by Holder's inequality from the bound 



l|S>.a(f)|| 



L 2 



< 



□ 



Now we prove part (a) of the proposition. Three tubes P,Q,R con- 
tribute to the integral only if they intersect. We consider the intersec- 
tion pattern of £ and £ + 77 tubes. For any £ tube Q, all £ + 77 tubes 
intersecting it are contained in a larger slab obtained by horizontally 
translating Q in the 77 direction 

H = 2(Q + {0}x [-277,217]) 
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We denote by Tt a locally finite covering of [— 1, 1] x R n with such slabs. 
Then we use the lemma to bound the integral in (58) by 

(PnQnR^ 

Using Cauchy-Schwartz in the second sum with respect to (Q, R) we 
bound this by 




/x5 (A + 1^.^)3 
where 

N = max \ {(Q, R); Q,RcH, P n Q n R ^ 0}| 

H,P 

An additional Cauchy-Schwartz allows us to estimate the above sum 
by 



N2 



'eimiO (EMM (ElMft 

\PeT v J \ q 



R 



^(A+|e-77|)2 

To conclude the proof it remains to establish a bound for N, namely 

N < m 

~ |{| + l«Afl| 

Both P C\ R and P C\Q have a A -1 time length and are uniquely deter- 
mined by this intersection up to finite multiplicity. It follows that 

N<\\I\ 

where I is the time interval where H and P intersect. Given the defi- 
nition of H it follows that / has the form 

/ = {t; \x + t£ + sr] | < 2 for some s G [—1, 1]} 

Symmetrizing we can assume that Xq = 0. Taking inner and wedge 
products with 77 it follows that t must satisfy 

1*1 < \t\ < IvW^vl' 1 

which lead to the desired bound for N. 

b) We first note that both Q and R are contained in spatial strips 
of size /j, x A oriented in the £ direction. Hence by orthogonality it 
suffices to prove the estimate in a single such strip. Then we can take 
advantage of the I 1 summability in the VX^ norm to further reduce 
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the estimate to the case when the r] tubes are spatially concentrated in 
a single /i x /i cube Z. 

A £ or a £ + rj tube needs a time of /iA -1 to move through such a 
cube. On the other hand, the tubes Q and R need a larger time 
to separate. Hence within Z we can identify the £ tubes and the £ + rj 
tubes. By orthogonality it suffices to consider a single Q and a single 
R. Consequently, the conclusion follows from the following counterpart 
of Lemma 5. 

Lemma 6. For £ and i] as in part (b) of the proposition the following 

estimate holds 

(61) 



Pcz 



^ u P v Q w R dx 

P€T V 



< V lluplli 



2^ \\up\\x v 



A2/Z 2 (A + |£ • 7]\)2 



The proof is almost identical with the proof of Lemma 5, the only 
difference is that we now have \I\ ~ fiX^ 1 . 

Case 2: Here we consider the product u^v\W\ where 1 < < A 
and prove a stronger bound, namely 



(62) 



u^vxwxdxdt 



< A ^KIUJMUJKIU* 



We use the modulation localization operators S^a/ioo to split each 
of the factors in two, 

Uf, = M <A 2 /100 m m + (1 - M <A 2 /100 M M ) 

etc. We observe that 



M^/wou^M^/woVxM^/woWxdxdt = 

due to the time frequency localizations. Precisely, the first factor is 
frequency localized in the region {\t\ < A 2 /50} while the other two 
are frequency localized in the region {r > A 2 /8}. Hence it remains to 
consider the case when at least one factor has high modulation. For 
that factor we have a favourable L 2 bound as in Proposition 11, 

and similarly for the other factors. Then it remains to prove that 

respectively 

WS^vxwx)]]^ < ^^IKIIxJhAlU;, 
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For this we need the Strichartz estimates in Proposition 10 (see also 
(32)). 

For the first bound we use the L « estimate for v\, respectively 
the L n+2 estimate for w M . 

For the second we first oberve that by orthogonality it suffices to 
prove it when both v\ and w\ are frequency localized to cubes of size 

//. Then we use the L ( « " estimate for both factors to derive L 4 
bounds by Sobolev embeddings. 

Case 3: Here we consider the product u^v\W\ where 1 < fi < A. 
This is treated exactly as above, and an estimate similar to (62) is 
obtained. 

□ 
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